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ABSTRACT 

The latest cosmological A^-body simulations find two intriguing properties for dark 
matter haloes: (1) their radial density profile, p, is better fit by a form that flattens to 
a constant at the halo center (the Einasto profile) than the widely-used NFW form; 
(2) the radial profile of the pseudo-phase-space density, p/of, on the other hand, 
continues to be well fit by a power law, as seen in earlier lower-resolution simulations. 
In this paper we use the Jeans equation to argue that (1) and (2) cannot both be 
true at all radii. We examine the implied radial dependence of p/of over 12 orders 
of magnitude in radius by solving the Jeans equation for a broad range of input p 
and velocity anisotropy j3. Independent of /3, we find that p/crj? is approximately a 
power law only over the limited range of halo radius resolvable by current simulations 
(down to ~ 0.1% of the virial radius), and p/crj? deviates significantly from a power- 
law below this scale for both the Einasto and NFW p. The same conclusion also 
applies to a more general density- velocity relation p/o~ e D . Conversely, when we enforce 
p/(Jr °c r 17 as an input, none of the physically allowed p (occurring for the narrow 
range 1.8 < rj ^ 1.9444) follows the Einasto form. We expect the next-generation 
simulations with better spatial resolution to settle the debate: either the Einasto profile 
will continue to hold and p/crj? will deviate from a power law, or p/of will continue as 
a power law and p will deviate from its current parametcrizations. 



1 INTRODUCTION 



iV-body simulations of cosmological structure formation 
have shown that the spherically-averaged radial profiles of 
the mass density and velocity dispersion of dark matter 
haloes follow simple and nearly universal functional forms 
that are largely independent of halo properties such as mass, 
environment, and formation history. For the density profile 
p, there had been considerable discussion about the value of 
the logarithmic slope of its central cusp, 7 = dlogp/dlogr, 
whether it is —1 as in the forms of, e.g.. 
|Navarro et~aL"1 ( |1997[ ), or 



McLa ughlin 20 05) |Hoffman et al.|[2007| |Stadel et aLl|2008l 
Navarro et al.|2008 1 . Both the total velocity dispersion and 
the velocity dispersion in the radial direction have been used 
to define a. Unlike the controversial density profile, whose 
best-fit form has changed over the years with improved nu- 
merical resolution, the power-law profile of p/cr has with- 
stood the scrutiny, and different studies have all reported 
similar findings except for a minor variation in the actual 
value of the slope of the power law. 



Hernquist (19901; 



Moore et al. (19991. Re 



-1.5 as in 

suits from recent TV-body simulations now suggest a lack of 
a definite inner slope - the density profile of these better- 
resolved dark matter haloes continues to flatten with shrink- 
ing radius (e.g.. 



A third relation 
anisotropy, (3(r) = 1 — at /a 
into account 



was proposed 

2 I J2 



when the velocity 
of simulated haloes was taken 



Hansen & Moore (2006) advocated a linear re- 



Navarro et al. 2004 



Merritt et al. 2005 



2006}|Graham et a!||2006| |Navarro et al.|2008| |Stade et al 



lation between /3 and the local logarithmic slope of p. Other 
studies, however, have found a large scatter in this relation, 
particularly in the outer parts of the haloes beyond the scale 



2008 1 . Functional forms such as Einasto ( 1969 1 and Prugniel 



6 Simien ( 1997 1 motivated by the Sersic profile for the sur- 
face brightness of galaxies ( Sersic|1968 1 appear to provide a 
more accurate fit to the latest simulations. In these forms, 

7 itself is a function of radius and asymptotes to zero at the 
halo center. 

A second property that has attracted much attention 
lately is the radial profile of the pseudo-phase-space den- 
sity, p/cr 3 , which has been reported to be well approximated 
by a power-law in a number of iV-body simulations (e.g., 
|Taylor fc Navarro|[200T| |Ascasibar etldT1|2"004l |Dehnen fc| 



radius ( |Dehnen fc McLaughlin|2005| |Navarro et al.|2 008 ) 

To help elucidate the physical meanings of these em- 
pirical relations determined from simulated haloes, a typ- 
ical approach is to use the Jeans equation for a spherical, 
self-gravitating collisionless system in equilibrium to predict 
the density or velocity structures of dark matter haloes un- 
der a certain set of assumptions. Most of the recent studies 
based on this approach have begun with the assumption of a 
power-law p/<r 3 (r). These papers then explored the density 
profiles allowed by the Jeans equation with either isotropic 
velocities (e.g. |Taylor fc Navarro|200ij|Hansen|2004||, or an 
anisotropic velocity profile j3(r) (e.g., Dehnen & McLaughlin 



© 0000 RAS 



2 Ma, Chang, & Zhang 



2005 Hansen Stadel 2006). Alternatively, some authors have 



chosen to study the radial profile of the velocity anisotropy 
P(r) from the Jeans equation, starti ng with a powe r-law 
p/u't and some input form for p (e.g., Zait et al.|2008 |. 

In this paper, we reexamine this question by taking a 
complementary approach. In Sec. 2, we do not assume a 
power-law p/a^> but instead solve for p/of using the Jeans 
equation with input p(r) and /3(r) that are motivated by 
simulations. We find that while p/a't is approximately a 
power law within the radial range probed by the current TV- 
body simulations, it is not universally a power-law and shows 
significant deviations below this range. This approach (for 
isotropic velocities) was also used in |Graham et aL"| f2006). 
While our results for the /3 = case agree with theirs in the 
range of radii explored in their paper (down to ~ 0.1% of 
the virial radius r v i r ), we explore 7 orders of magnitude be- 
low this range and predict non-power-law behavior in p/af, 
in the upcoming simulations. As we show in Sec. 3, this con- 
clusion remains unchanged when velocity anisotropy is in- 
cluded in the calculation, and when a more general form of 
the density- velocity relation, p/afj, is considered. For com- 
pleteness, Sec. 4 compares the results for the case when p/aj? 
is restricted to be a power law. 



a power-law in r: din p/dlnr — — 2(r/r_2) a . Unlike equa- 
tion (J2J that has a definite inner slope of 7, this profile con- 
tinues to flatten towards the halo center. Equation |3| has 
the same form as the Sersic profile commonly used to fit 
the two-dimensional projected surface brightness of galaxies 
l Sersic|1968) . The Sersic index n is simply equal to 1/ct. 

Starting with either profile in equation |2| or (|3|, and 
P = or some form of /3(r), we integrate the Jeans equa- 
tion |l]) to obtain oy(r). The numerical integration is per- 
formed from large radius (typically r = 10 5 r_2) down to the 
radius r of interest, using a standard Runge-Kutta integra- 



tor (Press et al. 19921. The velocity dispersion is assumed 



to be zero at the starting large radius. The oy from the in- 
tegration is then combined with the input p(r) to obtain 
p/al{r). 

The left two panels of Fig. 1 shows the result p/aj? (up- 
per panel) and its logarithmic slope (lower panel) over 12 
orders of magnitude in halo radius for (3 = and both types 
of input density profiles: GNFW with 7 = 0.5, 0.75, 1.0, and 
1.5 (dashed curves), and Einasto with a = 0.12,0.16, and 
0.18 (solid curves). The range of a is chosen to span the 
best-fit values of 0.115 to 0.179 for the six simulated galaxy- 



size haloes in Navarro et al. (20081. The right two panels of 



2 RESULT: p/a'f IS NOT A POWER-LAW 

For a spherical, self-gravitating collisionless system in equi- 
librium, its density and velocity structures obey the Jeans 
equation 



1 {pa 2 r ) i 2al(5 
P 



dr 



+ 



d$ 
dr 



(1) 



Fig. 1 show zoom-in views of the portion of the left figures 
that is resolvable by current simulations: 0.01 < r/r_2 < 10. 

Fig. 1 illustrates that p/a% is not a power-law in r for 
any of the seven input density profiles. For GNFW haloes, 
the slopes of p/aj? exhibit oscillations and deviate notice- 
ably from the critical case p/af. oc r -1 ' 9 (indicated by light 
dotted straight lines), in particular in the extreme cases of 
inner cusps of 7 = 0.5 (red dashed) and 1.5 (blue dashed). 
The Einasto haloes also deviate strongly from a power-law 
at small radius. The zoom-in panels show, however, that 
the slopes of p/af happen to be quite close to —1.9 over the 
limited range of r/r_2 ~ 0.01 to 10 that is resolvable by 
current simulations. This feature is particularly striking for 
the Einasto profiles, where all three solid curves for p/af, 
have a slope within —2.0 and —1.8 from r/r_2 ~ 0.01 to 10, 
with the deviations only starting to show up at the small- 
est radius r/r_2 ~ 0.01 near the simulation resolution limit. 
It is therefore not surprising that the power-law behavior 
of p/af(r) continues to appear to be valid even though the 
latest simulations find the Einasto form a better fit for p(r) 
than GNFW - Fig. 1 shows the Einasto profiles in fact pre- 
dict a more power-law p/af(r) for r/r_2 > 0.01 

The important point to note, however, is at the smaller 
radius of r/r_2 < 0.01 in Fig. 1. Here p/af(r) deviates far 
away from a pure power law with a wide range of slopes that 
depend on the input p. For Einasto haloes, the shape of p/af 
flattens continuously towards the halo center, reaching the 
(—1.5, —3) by |Moore et al.| ( |1999| . We find very similar re- asymptotic value of d\n p/af /d\nr = at r = regardless 



where p is the radial density profile, f3 — 1 — of /ay 3 is 
the velocity anisotropy parameter, ay and at are the one- 
dimensional radial and tangential velocity dispersions, and 
$ is the gravitational potential. To apply the Jeans equa- 
tion, we take as input a form for p and solve for ay. For the 
purposes of this paper, we consider two broad types of radial 
density profiles. The first type is cuspy all the way to the 
halo center with an inner logarithmic slope 7 and an outer 
slope 700: 



p(r) 



po 



(r/r- 2 )~>[l + r/r-2}' 



(2) 



where r_2, often referred to as the scale radius, is the 
radius at which din p/dlnr = —2. Examples of special 
cases of (7,700) that have been proposed for dark matter 
haloes include (-1, -3) by|Navarro et al.| (|1997[), (-1,-4) 



by Hernquist ( 1990 1 and Dubinski & Carlberg 



19911, and 



suits from our calculations for 700 = —3 vs —4; we will thus 
set 700 = —3 and refer to equation ([2| as GNFW below. 

The other type of density profile considered in this pa- 
per has a non-cuspy inner profile, given by the Einasto pro- 
file (Einasto|1969[) advocated in several recent studies (e.g., 



Merritt et al.|2005[|Gr"aham et al.|2006l|Navarro et al.|2008[ 



Stadel et al.|2008 1 



of the parameter a. This is not unexpected of the Einasto 
profile as the density approaches an asymptotic value in the 
core. For GNFW, there are two possibilities. For inner slopes 
of p that are steeper (shallower) than the critical 7 « 0.75, 
the power-law slopes of p/af are steeper (shallower) than 
the critical -1.9. At the critical 7 « 0.75, the GNFW halo 



has p/af 



at small r, a result consistent with that of 



ln^ = -[l-(r/r_ 2 ) c 

P-2 Ol 



(3) 



This profile has the feature that its logarithmic slope is itself 



Taylor & Navarro (20011, which showed that starting with 
an exact power-law of p/af oc r -1 ' 875 , the resulting density 
profile has an inner slope of dlnp/dlnr « 0.75. It is worth 
noting, however, that even the 7 = 0.75 GNFW halo shows 
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Figure 1. Radial profiles of the pseudo-phase-space density p/a^(r) (upper panels) and the corresponding logarithmic slope 
d In p/a^(r)/dlnr (lower panels) obtained from the spherical Jeans equation with = for seven input halo density profiles: Einasto 
(solid) with a = 0.18 (blue), 0.16 (green), and 0.12 (red), and GNFW (dashed) with 7 = 1.5 (blue), 1 (black), 0.75 (green), and 0.5 
(red). The left panels show the behavior of p/a^(r) over 12 orders of magnitude in r, while the right panels show zoom-in views of the 
region 0.01 < r/r—2 < 10, which corresponds to the range resolvable by the latest TV-body simulations. For ease of comparison with 
a power-law, the light dotted straight lines indicate the critical case p/crj? oc r~ 19 , and the y-axis in the upper right panel plots the 
logarithm of the ratio of p/a'f(r) to p/crf oc r~ 19 . All curves are scaled to have p/crj? = 1 at r = r—2- 



wiggles in the corresponding pja'l profile; that is, no GNFW 
haloes have an exact power-law p/o'f. 



3 FURTHER CONSIDERATIONS: p/c* IS 
STILL NOT A POWER-LAW 

We also find the conclusion reached in Sec. 2 to hold not just 
for isotropic velocity distributions, but also for anisotropic 
velocity distributions. To illustrate this point, we solve equa- 



tion (1) using an input /3(r) motivated by JV-body simula- 
tions ( jHansen fc Moore|2006] |Zait et al.|2008[ ), where (3 is a 
function of the local logarithmic slope of the density profile: 

We then compute p/a^(r) for a similar suite of GNFW and 
Einasto profiles. The results are shown in Fig. 2. What is 
especially notable is how insensitive the slopes of p/a 3 are to 
the form of /3(r) used in the calculation. This independence 
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-8 -6 -4-2 2 
Iog 10 r/r_ 2 



Figure 2. Effects of velocity anisotropy on the radial profile of 
p/erj? computed from the Jeans equation: (3 = (solid for Einasto; 
long dashed for GNFW; same as lower left panel in Fig. 1), and 
f} given by eq. Q (short dashed for Einasto; dashed-dotted for 
GNFW). The seven input density profiles are the same as in 
Fig. 1. This figure illustrates that including velocity anisotropy 
in the Jeans equation does not change p/cr'f significantly and still 
results in a non-power-law p/of. 




I°gl0 r A-2 

Figure 3. Radial dependence of three examples of the general 
density-velocity relation p/o"Jj (Sc hmidt et al,||2008[ l obtained 
from our calculations: p/cr 3,15 (solid), p/a^° (short-dashed), and 
p/of 18 (long-dashed). The y-axis shows the logarithmic slopes of 
these relations. The input (3 is from equation (4), and the input p 
are the three Einasto profiles shown in Figs. 1 and 2. This figure 
illustrates that like p/ff 3 , p/a^ is also far from being a power-law 
over the wide range of radius shown. 



from velocity anisotropy strengthens our conclusion reached 
earlier, namely, that p/a 3 is not a power law. 

Our findings are also in line with some recent work that 
has called into question the universality of p/o" 3 . For in- 
stance, |Schmidt et al.| ( |2008| > has advocated that individual 
simulated haloes are better fit by a generalized power-law 
relation that is not necessarily p/cr 3 : 

4cxr- a , (5) 



where an = oyVl + Dp, and D parameterizes a generalized 
an', for instance, D = 0, —2/3 and —1 correspond to o\d = 
(1-d), and a t , respectively. Schmidt et al. (20081 



, O"tot 



showed that the best-fit values of (D,e, a) differ from halo 
to halo, and as a set, they roughly follow the linear relations 
e = 0.97 D + 3.15 and a = 0.19D + 1.94. For a = oy (i.e. 
D — 0), the optimal relation is p/<r 3 15 oc r -1 ' 94 , which is 
consistent with the reported behavior of p/cr 3 m TV-body 
simulations within error bars. However, few haloes' best-fit 
value of D in Schmidt et al. (2008) is near D = 0. 



To assess whether any of these relations is closer to 
a power-law than p/o" 3 in our calculations, we choose 
three sets of D and e from their linear relation, (D, e) = 
(0, 3.15), (-2/3, 2.50), and (-1, 2.18), and plot in Fig. 3 the 
logarithmic slopes of these three relations p/crfj, using our 
solutions of the Jeans equation with non-zero j3 shown in 
Fig. 2. For clarity, only the three Einasto profiles are shown 
in Fig. 3, although our conclusions apply to the GNFW pro- 
files as well. The three solid curves in Fig. 3 (for D = 0) rep- 
resent p/cr 3 ' 15 , and are therefore very similar to the short- 
dashed curves for the Einasto p/of in Fig. 2. The other two 



sets of curves in Fig. 3 suggest that p/ofot (i.e. D — —2/3) 
and p/of' 18 (i.e. D = —1) are also far from being a power- 
law over the wide range of radius shown. When it is limited 
to the range of 0.01 < r/r_2 < 10 probed by simulations, 
p/ofot appears to be slightly closer to a power-law than p/a^ 
for the a = 0.12 (red) and 0.16 (green) Einasto profiles. Over 
the larger range of radius shown in Fig. 3, however, our ear- 
lier conclusion of a non-power-law p/o-j? carries over to the 
broader range of p/crf, shown. 



4 FORCING p/cr, 3 TO BE A POWER LAW 

Thus far we have solved the Jeans equation assuming an 
input p(r). As a point of comparison, we have also solved 
the Jeans equation following the works of previous authors 
with a starting assumption of a power law p/a'f. oc 



Most notably, |Taylor fc Navarro (20011 presented results 
for the special case of rj = 15/8 = 1.875, whereas Dehnen 



& McLaughlin ( |2005| ) (see also Austin et al. 2005 1 argued 
that only a single realistic solution exists for p, given when 
rj takes the particular value of 35/18 = 1.9444 for isotropic 
velocities, and r\ = 35/18 — 2/3o/9 for anisotropic velocities, 
where /3 = 0(r = 0). 

Here we examine a range of r\ and show in Fig. 4 our 
numerical solutions (for [3 = 0) from the Jeans equation 
for p(r) for four values of input power-law p/erf. oc r~ v : 
r\ = 1.8,1.875,1.9444 and 2.0. As r/ moves away from the 
critical value 1.9444, p drops to zero rapidly at some small 
radius if rj > 1.9444 (short-dashed curve), while p is cut- 
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Figure 4. Solutions for p(r) computed from the spherical Jeans 
equation with an input power-law p/af. oc r~ 71 and = 0. Four 
values of r\ are shown. Only when rj is in the range between ~ 1.8 
and 1.9444 does the solution not contain unphysical holes at small 
r (e.g. rj = 2) or excessive cutoffs at large r. The shapes of p in 
this narrow range of ??, however, deviate significantly from both 
the Einasto (solid red) and 7 = 1 GNFW (solid blue) profiles 
at small radius. Dark matter haloes therefore cannot follow the 
Einasto p and a power-law p/af simultaneously at all radii. 



off sharply at some large r if rj < 1.9444 (long-dashed and 
dotted curves). The former solution with a central hole is 
unphysical for dark matter haloes. The latter, however, is 
not automatically ruled out. Only for rj < 1.8 do we find 
the outer p to drop off too steeply to represent realistic 
ACDM haloes. It therefore appears that the narrow range 
of 1.8 < 77 ^ 1.9444 may admit physical solutions for p. 
It is important to keep in mind, however, that p in these 
cases are not well described at small radius by either the 
Einasto profile advocated by recent simulations, or the orig- 
inal 7=1 GNFW profile. Lowering the inner slope of the 
GNFW profile to 7 w 0.75 provides a better fit. Dark matter 
haloes therefore cannot be well fit by the Einasto p and a 
power-law p/af simultaneously at all radii. 



5 CONCLUSIONS 

Motivated by the apparent power-law radial profile of the 
pseudo-phase space density, p/af, reported in various N- 



body simulations (e.g., 


Taylor & Navarro 


2001 


Ascasibar 


ct al. 2004; Dchnen & McLaughlin 2005; Hoffman ct al. 2007; 


Stadel et al.||2008 


Navarro et al.||2008|), we solve the Jeans 



equation to study the implied pseudo-phase-space density 
for given parameterizations of p suggested by A r -body simu- 
lations, i.e., Einasto and GNFW. We find that independent 
of the velocity anisotropy, p/af is n °t a pure power law in 
radius for either the Einasto or GNFW profiles (left panels 
of Fig. 1). In the radial ranges that are probed by current 



A-body simulations (down to ~ 10 _3 r„i r ~ 10 _2 r_2), we 
find that p/af happens to be approximately a power law, 
in particular for the Einasto profile (right panels of Fig. 1). 
For radial scales right below the resolution of current simu- 
lations, however, we see significant deviations from a power 
law profile for p/a'r if either Einasto or GNFW continues 
to hold as a suitable parameterization of p. This result is 
unchanged when velocity anisotropy is included in the cal- 
culation (Fig. 2), and when a more general density- velocity 
relation of p/a^ is considered (Fig. 3). The hope to gain 
deep insight into the process of dark matter halo formation 
using a simple and universal power law scaling of p/af may 
therefore be misleading. Conversely, when p/af. oc r - '' is 
assumed as an input in the Jeans equation, none of the real- 
istic solutions for p, which occur only for the narrow range 
of 1.8 < 77 < 1.9444, take the Einasto form (Fig. 4). 

We therefore conclude that the two halo properties seen 
in current simulations - the Einasto p and power-law p/af, 
— cannot hold simultaneously at all radii. We expect that 
the upcoming simulations with better spatial resolution will 
settle this debate: either the Einasto profile will continue to 
hold and p/af, will show a break from a power law, or p/af, 
will continue as a power law inward and p will deviate from 
its current parameterizations. 

We thank Onsi Fakhouri, Peng Oh, and Simon White 
for useful discussions. PC and JZ are supported by the The- 
oretical Astrophysics Center at UC Berkeley. 
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